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ABSTRACT 

The role of aliphatic side groups on the formation of astronomical uniden¬ 
tified infrared emission (UIE) features is investigated by applying the density 
functional theory (DFT) to a series of molecules with mixed aliphatic-aromatic 
structures. The effects of introducing various aliphatic groups to a fixed poly¬ 
cyclic aromatic hydrocarbon (PAH) core (ovalene) are studied. Simulated spectra 
for each molecule are produced by applying a Drude profile at T=500 K while 
the molecule is kept at its electronic ground state. The vibrational normal modes 
are classified using a semi-quantitative method. This allows us to separate the 
aromatic and aliphatic vibrations and therefore provide clues to what types of 
vibrations are responsible for the emissions bands at different wavelengths. We 
find that many of the UIE bands are not pure aromatic vibrational bands but 
may represent coupled vibrational modes. The effects of aliphatic groups on the 
formation of the 8 /irn plateau are quantitatively determined. The vibrational 
motions of methyl (—CH 3 ) and methylene (—CH 2 —) groups can cause the merg¬ 
ing of the vibrational bands of the parent PAH and the forming of broad features. 
These results suggest that aliphatic structures can play an important role in the 
UIE phenomenon. 

Subject headings: ISM:lines and bands — ISM: molecules — ISM: planetary 
nebulae: general — infrared: ISM 
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Introduction 


While it is widely known that interstellar dust contains an aromatic component, the 
presence of an aliphatic component is less well appreciated. The 3.4 /irn feature was f irst de¬ 


tecte d in absorption against the infrared background of the galactic center (jSoifer et al. 


1976h. but_ no identificat i on wa s made at that time. This feature was rediscovered by 


Wic kramas in ghe fe Allenl (11980 1 a few years later in a search for organic dust in the dif¬ 


fuse interstellar medium. Higher spectral resolution observations of the 3.4 /irn feature in 
the planetary nebula IRAS 21 282+5050 show that it is compo sed of several subfeatures at 


3.4, 3.46, 3.51, and 3.56 /zm (iJourdain de Muizon et al.l Il990h . which can be identified as 
C—H stretching modes of methyl (—CH?) and methylene (—CHA groups in aliphatic hydro- 


carbon materials ( 

Dulev & Williams 

1983) 

. Since then, the feature has been seen in emission 

in proto-i 
galaxies ( 

planetary nebulae 

Geballe et al. 

19921; 

Hrivnak et al. 20071 and in absorption in 

Wright et al. [l996; 

Dartois et al. 

2004, 

2007 

)• 


In addition to the 3.4 /irn aliphatic C—H stretching mode, features at 6.9 and 7.3 /irn 
corresponding to the bending m odes of aliphatic methyl and methylene groups have been seen 
as early as 1979 (IWillner et al.lll9791 ). Th ese features have also been de tected in absorption 
flChiar et al.l 120021) as well as in emission (IK wok. Volk fe Hrivnaklll999l) . 


Akari observations of M82 have shown that the strength of the 3.4 /irn feature increases 
from the disk to the halo of the galaxy, sugges ting that aliphatic organics are distributed 
throughout the volume of galaxies (iYamagishi et al. 2012). Since abundance can easily be 
derived from the strengths of absorption features, it is estimated that the amount of carbo n 
tied up in aliphatic compounds represent 15 — 30% of the total interstellar carbon (Dartois 

201lh . 


Another manifestation of the aliphatic component is the strong, broad emission plateau 
features seen in planetary nebulae and proto-planetary nebulae. The plateau features cen¬ 
tered around 8 and 12 /zm were first s een in the JR AS LRS (IKwok. Volk fe Hrivnakl 119891 1 
and the Kuiper Airborne Observatory ( Bu ss et al. 19901) spectra of proto-planetary nebulae. 
These two plateau features are much better defined by ISO observations and have been 
identified as su perpositions of a variety of aliphatic side groups attached to aromatic rings 
(Kwok, Volk & Bemathl l200ll ). These broad features cannot be explained by the traditional 
pure aromatic models such as the polycyclic aromatic hydrocarbon (PAH) molecules. The 
presence of these features, as well as the 3.4, 6.9, and 7.3 /irn features, suggests that the 
chemical structures of organic dust in the interstellar medium is much more complex than 
previously believed. 


By introducing H into graphite ( sp 2 ) and diamond (sp 3 ), a variety of amorphous C—H 








































































3 


alloys can be created (Rob e rtson 2002). Different geometric structures with long- and short- 
range can be created by varying the aromatic to aliphatic ratios. The spectral charac¬ 
teristics of these structures and their possible role as carrier of the u nident ified infrared 
(UIE) bands have recently been extensively studi ed (Jones! 


emission 


2012al 


315 ) • 


Labo¬ 


ratory studies of such amorphous c arbonaceous solids fjDischler. Bubenzer fe Koidll Il983al: 
iDischler. Bubenzer. fe Koidll Il983bh suggest that they show many similar spectral charac¬ 
teristics as astronomical UIE bands. 


The vibrational pro perties of PAH molecules have been extensively studied fe.g-. lBanschlicher fe Langh 
1997|; Ricca et al. 120131) . Infrared spectra of PAHs with different geometries, sizes, and 


charge states have been obtained with quantu m chemical calculations flMalloci et al.l 120071: 
Bauschlicher et al.l 12010 : Boersma et ah 2014). Recently, there has been an increasing in¬ 
terest in the effects of aliphatic structures. The UIE bands have been sug gested to be 


better explained by mixed aromatic/aliphatic organic nanoparticles (MAON, Kwok fe Zhang 


2 0111) . The contributions of methyl and methylene groups have been computed by Yan g et al. 


(2013), and the vibration al s_pectra of a variety of aliphatic structures have been calculated 
by Papoular & Papo ul ar (2014). 


Through advances in computational chemistry methods and the increasing power of 
computing hardware, it is now possible to perform theoretical investigation of the vibra¬ 
tional properties of large organic molecules (Papoular 12013 : Ricca et al. 120131 ). With such 
theoretical studies, we can identify the explicit roles of aliphatic groups in the vibrational 
spectrum of the molecules. In this work, we aim to study the changes in the infrared emission 
spectra of a PAH molecule upon addition of the aliphatic groups and to specifically identify 
the types of vibrational motions responsible for such changes. 


2. Computational methods 


Our calculations are b a sed on the De nsity Functional Theory fDFT. lHohenberg fc Kohn 


1964: iKohn fc Shaml Il965l: IZieglerl Il99ll) where the molecular non-relativistic Schrodinge r 
equation is solved under the Born-O ppenheimer approximation flBorn fe Oppenheimerlll927l ). 
The BH&HLYP hybrid functional (iBeckelIl993t ) in combination with polarization consistent 
basis set PCI (jJensen 200ll , 20021) are applied to obtain the equilibrium geometries and 
the fundamental vibrational frequencies of hydrocarbon molecules. The chosen PCI basis 
set includes the necessary polarized d, and p orbitals on C and H centers individually. The 
accuracy of the vibrational f requencies under the BH&HLYP/PC1 model is estimated to 
have an rms error of 2 cm -1 OLaurv et al.l 2012). All geometries have been optimized under 
the default criteria, i.e. forces < 4.5 x 10 -4 hartree/bohr and rms error of the forces less than 
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3.0 x 10~ 4 . Fine integration grid in combination with the default convergence criterion of 
1 x 10” 8 is applied for self-consistent field calculation of DFT density matrix. The optimized 
geometries are all characterized as local minima, established by the positive values of all 
frequencies and their associated eigenvalues of the second derivative matrix. The_vibrational 
frequencies are calculated via the harmonic normal mode approximation (iOchterskil 19991 ). 
The scaling factors of 0.9311 and 0.9352 fop vibrational frequencies > 1000 cm -1 and < 
1000 cm -1 , respectively, are then applied flLaurv et al. 2012). The final infrared spectra for 
each isomer are constructed by plotting the scaled frequency values versus the calculated 
vibrational mode intensities. 


The simple quantitative decomposition normal mode analysis is adopted in combination 
with customary visualizing of the normal modes. The Cartes ian displacements x, y, z of 
all the atoms calculated in the vibrational analysis flOchtersk i 119991 ) are converted to a 
quantitative measure R , which is defined as 


R = \Jx 2 + y 2 + z 2 (1) 

where x , y , z, R are all expressed in units of A. 

The analysis is done with a small code called Vibanalysis written by SAS. This code 
has the ability to read the atoms displacement in cartesian coordinates and to perform the 
described analysis. Its output is the list of the atoms (labeled as aliphatic or aromatic) with 
the major contribution in the corresponding normal mode, i.e. the atoms with R > 0.1 A. 


All ab i nitio c alculations have been performed by Gaussian 09, Revision C.01 package 
(Frisch et ah 2009). The HKU grid-point supercomputer facility equipped with 8 and 12 
CPU cores and 30—40 GB of RAM on each node has been utilized for this purpose. Visual¬ 
izing and manipulating the results of vibrational normal mode analysis were performed by 
utilizing the Chemcraft suit program 


In order to facilitate comparison with the observed astronomical infrared emission spec¬ 
tra, the calculated vibrational-mode intensities have been converted to a simulated observed 
spectra by applying a thermal excitation and a Drude profile. Assuming a thermal popula¬ 
tion, the intensities of emission lines can be calculated from the vibrational-mode intensities 
calculated by DFT. We assume that the intensity of the emission line at certain wavelength 
is equal to its absorption intensity calculated by DFT. Ignoring the contributions from hot 
bands and overtones, the line intensity of the ith vibrational mode (/,;) is proportional to 
exp{—hvi/kT) where zzj is the frequency of the normal mode, and h and k are the Planck 


2 http: //www. chemcraftprog.com 
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and Boltzmann constants, respectively. A Drude profile 


H A) = 


2/, A, 


x 


7 


( 2 ) 


7TC7 (A/A; - A;/A) 2 + y 2 
is applied to broaden all the IR features, where A; is the central wavelength, 7 is the fractional 
full-width at half maximum (FWHM) of the simulated peaks and c is the speed of light. The 
value of 7 chosen to fit the experimental or observational spectra. The modeled spectra are 
then obtained by sum over all F;(A) values: 


F{ A) = EF;( A) (3) 

where F( A) is the final intensity of the emission band at each point on the Drude profile. In 
this model no electronic excitations have been considered therefore all the molecules in our 
present study are kept in their singlet electronic ground state. T = 500 K is chosen for the 
simulation tempe r ature of emission spectra as this is the optimum temperature suggested by 
Cook fe Savkallvl (119980 that keeps the long wavelength vibrational levels well populated. 


3. Theory versus experiment 


I11 order to check the accuracy of our quantum chemical model in predicting the frequen¬ 
cies and intensities of the fundamental vibrational normal modes in hydrocarbon molecules, 
we compare our theoretical IR spectra with the re cent expe rimental FT-IR data on num¬ 
ber of hydrogenated PAH molecules reported by ISan dford et al.l (1201311 . These are the 
Group H molecules: C 2 2 H 16 (DHDB[ah]A), C 2 oHi 4 (DHB[e]P) and C 20 H 14 (TRIP), the largest 
molecules in the set. Qualitatively we found very good agreement between the calculated 
bands position and intensity ratios and their corresponding experimental counterparts in the 
IR spectra. A comparison between the calculated IR spectrum with the experimental spec¬ 
trum of TRIP is shown in Figure CD We have also evaluated the accuracy of our theoretical 
results quantitatively and the results are shown in Table 1. The calculated spectra are broad¬ 
ened using the Drude function with the FWHM equal to that of experimental spectra. Each 
peak is actually a composition of multiple transitions. The peak positions and integrated 
strengths are measured and listed in Table 1. We find that the BH&HLYP vibrational modes 
intensity ratios are in very good agreement with the reported experimental band strengths. 
Accuracy of less than 0.2 /irn is achieved for peak positions of the bands with the the largest 
theoretical error found for 8.8 /irn band in Table 1. 

The reliability of our method to simulate the infra red e mission spectrum was verified 
against experimental spectra of PAH molecules (Cook et al. 1996]). Comparisons are made 
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with three selected PAH molecules, phenanthrene (C 14 H 10 ), pyrene (CigHio) and coronene 
(C 24 H 12 ) (Figure [2]). While there are differences in relative intensities, the positions of 
the theoretical bands are in good agreement with experimental values. The average absolute 
error of the position of simulated bands are 15.31 cm" 1 , 20.68 cm " 1 and 20.47 cm ' 1 (equivalent 
to 0.12 to 0.13 /im error) for phenathrene, pyerene and coronene, respectively. The 10.42 /mi 
band in exp erimen tal spectra of coronene (Figure [2]) is a contamination band as reported by 


Cook et ah (119981 ). 


4. Interpretation of emission spectra 


The detection of astronomical aliphatic emission features has led to interest in the ex¬ 
ploration of variations from the PAH model, which until recently has been the most widely 
accepted model for the astronomical UIE bands. Recent efforts include the laboratory study 
of partially hydrogenated PAHs ( Sandford e.t al. 2013) and theoretical studies of the vibra¬ 


tional modes of kerogen derivatives ( Papoular 20131 ). In this paper, we explore the interme¬ 


diate case through the gradual addition of aliphatic side groups to an aromatic core. We 
begin with a medium-size PAH molecule, ovalene (C 32 H 14 , Figure [3]), and gradually replace 
its peripheral H atoms with the different —CH=CHCH 3 , —CH 2 CH=CH 2 , —CH 2 CH 2 CH 2 —, 
—CH 2 (CH 2 ) 3 CH 3 , —C(CH 3 ) 3 and —CH=C(CH 3 ) 2 aliphatic functional groups. Such struc¬ 
tures are similar to those described in Figure 4 of Kwok, Volk V Bcrnath (2001). With 
different type of functional groups the symmetry of the molecule is reduced to C\ and there¬ 
fore the number of IR active modes is maximized. The chemical formula, 2D structural 
formula, 3D ball and stick model for the local-minimum geometries of our sample molecules 
are shown in Figure [TJ 


In order to compare with astronomical spectra, we apply the Drude profiles to the 
molecules, creating broad emission bands from the vibrational lines. The normal mode 
vibrational line(s) that coincides with the peaks of emission band are used to classify and label 
the band. We End that in majority of the cases, the most intense IR lines under a particular 
emission band coincide well with the band peak position. For bands that overlap with similar 
strength intensity IR lines, we rely on emission peak position to pick the corresponding 
vibrational mode. 


The vibrational modes of the ovalene molecule in the 2—20 /mi range are identified in 
Table [2l Each of the vibrational mode is assigned based on our analysis on displacement 
vectors together with visual inspection of the animation of the vibrational modes. These 
results are generally con sistent with usual interpretation of vibrations for aromatic/aliphatic 
compounds (Duley & Williams 19831 ). For example, the 3.3 /mi band is assigned as aro- 
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matic C— H stretch. However, some other bands are the result of coupled motions and their 
vibrational mode assignments are not as straight forward as the 3.3 /irn feature. 

Figures I5HTU1 show the spectral changes as aliphatic components are added to the ova- 
lene molecule. A glance at these figures show that the PAH bands at 8.77 /irn, 11.78 /irn, 
13.04 /irn, 15.59 /irn and 17.94 /irn weaken relatively upon the addition of aliphatic groups. 
Qualitatively, all of these bands show the contribution from in-plane or out-of-plane bending 
vibrations of aromatic C— H bonds (Tabled]). The weakening of these bands is the result of 
replacing these hydrogens by aliphatic groups. 

Below, we discuss in more details the spectra in separate wavelength regions. 

• 3—4 /rm region (Figured]) 

The ovalene molecule has a single emission band at 3.25 /irn, which is due to aromatic 
C—H stretch. When aliphatic groups are added, a broader feature around 3.45 /irn 
begins to appear. Vibrational analysis shows that this band is due to aliphatic (sp 3 ) 
C—H stretching modes. The addition of more aliphatic groups causes the two bands 
to merge, as can be seen in the spectrum of C 55 H 52 . These results are in agreement 
with functional group frequency assignment sc heme widely used in organic chemistry 


(Fuhrer et ah 

1972S; 

Giinzlcr & Gremlich 

2002) 


• 5.6—7.6 /rm region (Figured]) 

The 6.2 /irn band is a common feature among all the molecules under study. Normal 
mode analysis shows that this band is due to aromatic C=C stretching coupled with in¬ 
plane bending mode of aromatic C—H bonds. An aliphatic feature is observed as a weak 
peak at 6.01 /irn at the shoulder to the aromatic 6.2 /irn band. This 6.01 /irn feature can 
be identified as aliphatic C=C stretching coupled with aliphatic =C-H in-plane bending 
modes found in H 2 C=CH—CH 2 — and H 3 C—CH=CH— functional groups. This band 
is blue shifted to 5.94 /irn when the aliphatic double bond is conjugated to aromatic 
skeleton (C 3 sH 22 ). In this case the carbon atom of the ring bearing the aliphatic group 
also possess the considerable displacement and thus our normal mode analysis counts 
this motion as one of the mixed modes (C 3 sH 22 ). Ovalene (PAH) shows a weak band 
at 6.59 /in 1 composed of coupled in plane C=C stretching and C-H bending vibrations 
(Table [2]). The band shifts to 6.64 /irn upon replacement of the first hydrogen atom by 
H 3 C-CH=CH- group (C 35 Hi 8 ). This band vanishes completely in all other molecules 
in our series. 

Aliphatic groups vibrational modes split the ovalene emission band at 7.19 /on into 
two bands. One centered at 7 /an and the other band at 7.25 /an. The vibrations at 








7 pm are found to be deformation motions of — CH3 (methyl) groups and scissoring 
motions of aliphatic — CH 2 (methylene) groups fjGiinzler fe Gremlichl 2002). These 
aliphatic motions combine with aromatic vibrational modes at this region to form 7.25 
/im band. Addition of more aliphatic —CH 3 and —CH 2 groups results in merging these 
two bands (7 pm and 7.25 pm) and the formation of the broadened band at 6.95 pm 
(C 55 H 52 ). This is also supported by the increase in the number of aliphatic and mixed 
modes as calculated via our normal mode decomposition analysis. 


7.6—9.6 pm region (Figure [7]) 

In this spectral region, we see a marked growth in the number of pure aliphatic and 
mixed modes. As aliphatic groups are added to the ovalene molecule, three distinct 
bands of ovalene at 7.80, 8.22, and 8.77 pm are merged together and appear as a one 
broad feature at around 8 pm in CssH^. Our vibr ational analy s is show s that the 
coupling of — CH 2 groups twisting vibrational modes fjGiinzler & Gremlichl 2002) with 
the aromatic vibrations at this region is the main reason for this change. Other types 
of aliphatic motions include H—C=C—H (para, in- plane bending), —CH 3 deformation, 
—C—H (aromatic, bending) modes. The strong 8.77 pm band in ovalene is blue shifted 
gradually upon combining with aliphatic vibrations and ultimately appears as weak 
band around 8.54 pm in the CssH^. 


The vibration under the weak 9.54 pm band of ovalene can not simply be assigned by 
usual terminology. We suggest using in-plane rocking mode of aromatic HC=CH unit 
for labeling this vibration. This band is replaced by the weak band in the range of 
9.1—9.3 pm as a result of combination with —CH 3 and —CH 2 deformation, wagging 
and twisting modes. The band intensity (9.1—9.3 pm) increases when these aliphatic 
groups are added in the form of cyclohexane group (C 47 H 36 ). 


9.4—11.4 pm region (Figured]) 

C 55 H 52 , the molecule with the largest number of aliphatic groups, shows a weak band 
at 9.7 pm . This band can be identified as coupling between various twisting/wagging 
modes of —CH 2 groups and deformation mode of — CH 3 groups and aromatic in-plane 
C—H bending modes. The vibration of aromatic C—H bonds are different from what 
we observed in the ovalene molecule at 9.54 pm, which weakens with the addition of 
aliphatic groups. 

The out-of-plane bending modes of para hydrogens in H 3 CHC=CH— group results in 
appearance of a band at 10.24 pm (see the spectrum of C 35 H 18 ). This feature exists for 
all the molecules in our series. One of the other characteristic vibrations that appear 
under the 10.24 pm band is the stretching mode of the C—C single bond in the pentane 
side chain. 
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The most notable feature in this spectral region is the strong band of ovalene (PAH) at 
11.11 pm (Table [2]). The combination of aromatic —C—H out-of-plane bending modes 
with deformation vibrations of — CH 3 groups shifts this strong band to 11.03 /jrn band, 
as can be seen in the spectrum of the C 35 H 18 molecule. Addition of H 2 C=CHCH 2 — 
aliphatic group causes the formation of another band at 10.85 pm. The aliphatic vibra¬ 
tions under this band are out-of-plane bending modes of olefin C—H bonds (H 2 C=). 
These vibrations are coupled with the out-of-plane bending modes of aromatic C—H 
bonds. By increasing the numbers of —CH 3 and —CH 2 groups, the deformation modes 
of — CH 3 groups and twisting or rocking modes of —CH 2 groups are added to the 
olehnic vibrations at this region of spectra. The coupling between all these aliphatic 
modes and the aromatic out-of-plane C—H bending modes results in the replacement 
of strong 11.11 pm feature in ovalene with the 10.92 /irn broad band in C 55 H 52 . 

• 11.4—14.0 /irn region (Figure [9]) 

The strong 11.78 pm band (Table [2]) in ovalene is gradually weaken by the decrease 
in the number of aromatic C—H bonds in the other molecules. At the same time, the 
band position is red shifted to 12 pm by combining the C—H out-of-plane bending mode 
with wagging and deformation modes of — CH 2 and — CH 3 , especially the methyl and 
methylene group on the cyclohexane ring. The coupling of various C—C—C aromahc bond 
angle deformation modes of aliphatic groups (for instance: H 2 C=CH—CH 2 -C aromat * c 
with the PAH modes at 13.04 pm, forms a very weak band at around 13.34 pm in 
C 55 H 52 . This feature can hardly be seen in Figure [U] 

• 14.0—20.0 pm region (Figure ITOl) 

While the bands of ovalene (PAH) can be identified by the usual vibrational modes 
(Figure [TT]), the bands of the other molecules in our series in this part of the spectrum 
cannot be so easily labeled (Figure [T^]). Our simple displacement vector analysis shows 
that the motions involve a large fraction of the molecules and are more complex in 
nature. However, distinction between aromatic and aliphatic contributions can still 
be made. The ovalene (PAH) molecule has a weak band at 14.45 pm and two strong 
bands at 15.59 pm and 17.94 pm. These bands are weaken by replacement of aromatic 
C—H bonds with aliphatic groups. The only noticeable fingerprint of the aliphatic 
group in this region is the concerted rocking motion of — CH 2 groups on the pentane 
side chain. This vibration forms a band at 15.31 pm in CssHs^ Based on our analysis 
it is suggested that the strong 17.94 pm band (Table [2]) in ovalene is blue shifted to 
16.9 pm region and then it gradually vanishes. 
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5. Discussions 


The UIE phenomenon is a complex one. In addition to the aromatic features at 3.3, 6.2, 
7.7, 8.6 and 11.3 fin i, there are aliphatic features at 3.4 and 6.9 fin i. There are also uniden¬ 
tified fea tures at 15.8, 16.4, 17.4, 17.8 , and 18.9 fim whic h are observed in p r o t o- p 1 a net a r y 
nebulae ( Kwok. Volk &; Hrivnakl 1999 ). reflection nebulae (jSellgren. Uchida fe Werner 2007 1. 


and galaxies (Sturm et ah 200011 . Broad emission plateau features at 8 and 12 /jrn, and 
maybe also around 17 fim are commonly observed. In recent years, there has been an 
increasing interest in the plateau features as they are found to b e associated with the 
presence of f ullerene fCsnl (IZhang fc KwokH2013t lOtsuka et al.l 1201411 . The unidentified 21 


fim feature (Kwok, Volk & Hrivnak 


19891 


is also often accompanied by the 8 and 12 fim 
plateau features flCerrigone et all 20 111: Volk et al. 2011). Suggestions have been made that 


the carrier of the plateau fea tures are precursors to fullerene (IBernard-Salas et al.l 12012 


Garcfa-Hernandez et ah 2012). A good understanding of the origin of the plateau features 


is key to the identification of chemical structure of the carrier of the UIE bands and the 
history of chemical synthesis in the late stages of stellar evolution. 


How well can a mixed aromatic/aliphatic model work in simulating the astronomical 
UIE bands in comparison to the PAH model? Traditiona 11 y.Alie UIE ba nds ar e explained by 
collective emissions from a mixture of PAH molecules. Cook fe Saykallv (119981) used a family 
of nine neutral PAHs (including ovalene) to simulate the astronomical spectra. They found 
that the inclusion of cation PAHs can explain the UIE better than a collection of neutral 
PAHs, but their simulation still cannot reproduce the details of UIE. They suggest that a 
larger number of PAHs and/or larger PAHs are required. With data o n hundr eds of PAH 
molecules now available from the NASA Ames PAH database (Boersma et ah 2014J), it is 
possible to use PAH molecules with different sizes and charged states to simulate the astro¬ 
nomical UIE spectra (iCamil 2011). However, questions have been raised on how meaningful 
such fitting exercises are due to the large number of parameters involved (Zhang & K wok 


2015). 


Although the range of mixed aromatic/aliphatic structures that we have calculated is 
limited, Figure [PlKC^H^. an ovalene molecule with aliphatic branches) shows that aliphatic 
side groups can also create a host of features in the 6—9 and 10—14 fim regions, and therefore 
offers the possibility of explaining the 8 and 12 /mi plateau emission features. In the 15—20 
fim region, the mixed aromatic/aliphatic structure can lead to complex vibrational motions, 
and could form the basis of explaining the 17 fim plateau emission feature. These results 
suggest that the inclusion of aliphatic side groups will help our interpretation of the UIE 
phenomenon. 





















































6. Conclusions 


We have conducted a theoretical study on the changes on the infrared emission spectra of 
the PAH molecule ovalene through the replacement of aromatic hydrogen atoms by aliphatic 
groups. The local-minimum geometry and the corresponding vibrational normal modes are 
calculated by applying the density functional theory. The IR absorption lines positions are 
then corrected for anharmonicity, couplings and other complex electronic effects, by applying 
the most recently developed vibrational scale factor scheme. In order make comparison with 
astronomical spectra, these absorption lines are translated to emission bands assuming a 
Drude profile and thermal excitation. The resulted emission spectra show good agreement 
with available laboratory experimental data. 

In order to understand the type of vibrational motions associated with each simulated 
emission bands, the calculated absorption lines are mapped back to the corresponding sim¬ 
ulated emission spectra. Through the use of a simple quantitative analysis of cartesian 
components of displacement vector of all atoms for each normal mode (IR or Raman active), 
we are able to separate the contributions of aromatic and aliphatic parts. Results show that 
the major regions of IR spectra affected by aliphatic groups are 3—4, 6—7, 7—8, and 8—9 
/jrn. Among the last three regions, the region of 7—8 /jrn contains the largest number of 
aliphatic vibrations. This is the region of spectra that 8 /jrn broad plateau feature is observed 
in astronomical sources and we suggest that aliphatic side groups play a major role in the 
formation of the 8 /irn plateau feature. 

The theoretical scheme described in this work can be very useful in the interpretation 
of the astronomical spectra and the identification of the structure of the molecular carrier of 
the UIE bands. In this paper, we have limited ourselves to the study of molecules with an 
aromatic core and aliphatic branches. In order to further explore the MAON model, we will 
study more irregular structures with random aromatic rings connected by randomly oriented 
aliphatic branches. Since such structures are more complex, larger computational resources 
are required. We hope this work will motivate further theoretical studies of large, amorphous 
organic compounds for a better understanding of the origin of the astronomical UIE bands. 
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Fig. 1. — The simulated IR spectra o f C 2 oH 14 (shown in blue) at 500 K compared with 
the experimental data of S andf ord et alj (120131 ) (shown in red). The broadening parameters 
were chosen to match to experimental broadening. The peak wavelengths of some of the 
features are also marked: theoretical values in normal font and experimental values in bold 
italics. The optimized geometry of the molecule is shown in the insert. (A color version of 
this figure is available in the online journal) 
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Fig. 2.— Comparison of the experimental spectra (in red) of phena nthrene (top panel), 
pyrene (middle panel), and coronene (bottom panel) flCook et al.l 19961) with our theoretical 
DFT spectrum (in blue). The numbers are peak wavelengths of the features with the ex¬ 
perimental values in bold italics and the theoretical value_s in regul ar font. The * marks a 
contamination band in the spectrum of coronene (Cook et al. 1998). The digitized experi¬ 
mental spectrum intensity values have been scaled up by a factor of 10 for ease of comparison. 
(A color version of this figure is available in the online journal) 
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Fig. 3.— Calculated 2—20 pm spectrum of C 32 H 14 . The 2D structure and 3D capped sticks 
representation of the optimized geometry of C 32 H 14 are also shown. 






Fig. 4.— The 2D and 3D capped stick representations of the optimized geometry of our 
sample molecules. The molecules are (from left to right and from top to bottom) C 35 H 18 , 
C 38 H 22 , C 42 H 28 , C 47 H 38 , C 51 H 46 , and C 55 H 52 . 
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Fig. 5.— Simulated 3.0—4.0 /jrn emission spectra of the 7 molecules with increasing degrees 
of aliphatic structures. The first molecule C 32 H 14 (ovalene) is a purely aromatic molecule.. 
(A color version of this figure is available in the online journal) 
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Fig. 6.— Simulated 5.6—7.6 pm emission spectra of the 7 molecules with increasing degrees 
of aliphatic structures. (A color version of this figure is available in the online journal) 
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Fig. 7.— Simulated 7.6—9.6 /irn emission spectra of the 7 molecules with increasing degrees 
of aliphatic structures. (A color version of this figure is available in the online journal) 
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Fig. 8.— Simulated 9.4—11.4 /irn emission spectra of the 7 molecules with increasing degrees 
of aliphatic structures. (A color version of this figure is available in the online journal) 
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Fig. 9.— Simulated 11.4-14.0 /im emission spectra of the 7 molecules with increasing degrees 
of aliphatic structures. (A color version of this figure is available in the online journal) 
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Fig. 10.— Simulated 14.0—20.0 [i m emission spectra of the 7 molecules with increasing 
degrees of aliphatic structures. (A color version of this figure is available in the online 
journal) 
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Fig. 11.— A snapshot of animation showing the vibrational motion associated with the 15.59 
/ini feature of ovalene (C 32 H 44 ). The direction and magnitude of motions are indicated by 
arrows. (A full video is available in the online journal) 
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Fig. 12.— A snapshot of animation showing the vibrational motion associated with the 15.53 
/irn feature of C 55 H 52 . The direction and magnitude of motions are indicated by arrows. (A 
full video is available in the online journal) 
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Fig. 13.— Comparison of the simulated emission spectrum of C 55 H 52 (solid line) with the 
Infrared Space Observatory spectrum (dotted line) of the planetary nebula BD+30°3639 (in 
units of F\). The narrow features in the astronomical spectrum are atomic lines. 
















Table 1. 


The comparison between simulated IR spectrum and experimental FT-IR data. 



Experiment® 


DFT 

Absolute error 6 
(/mi) 

Position 
(cm -1 ) (/<m) 

Strengths® 

Position 

(/mi) 

Strengths® 

626.2 

15.96 

0.57 

16.07 

0.35 

0.11 

742.4 

13.47 

1.0 

13.51 

0.70 

0.04 

796.7 

12.55 

0.15 

12.59 

0.15 

0.04 

1166.2 

8.57 

0.11 

8.76 

0.17 

0.19 

1193.2 

8.38 

0.19 

8.44 

0.24 

0.06 

1460.4 

6.85 

0.66 

6.90 

1.0 

0.05 

2958.3 

3.38 

0.2 

3.35 

0.13 

0.03 

3028.6 

3.30 

0.84 

3.25 

0.23 

0.05 


a The experimental data are the peak positions and strength of the dominant 
band at each region of spectra (Sandford et al. 2013). 


1A, 


exp 


A 


theo | 


Normalized to the strongest line. 
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Table 2. Identification of the vibrational modes of emission bands of ovalene. 


Bands (/mi) 


Vibrations 


3.26 

6.18 

6.59 
7.19 
7.80 
8.22 
8.77 
9.54 
11.11 
11.78 
13.04 
14.14 

15.59 
17.94 


C-H stretching 

C=C stretching coupled with in-plane C-H bending 
C=C stretching coupled with in-plane C-H bending 
C=C stretching coupled with in-plane C-H bending 
C=C stretching coupled with in-plane C-H bending 
C-H in-plane bending 
C-H in-plane bending 

[C=C] in-plane rocking mode unit coupled with in-plane C-H bending 
C-H out-of-plane bending coupled with out-of-plane ring deformation 

C-H out-of-plane bending 

C-H out-of-plane bending / In-plane ring deformation 
In-plane ring deformation 

C-H out-of-plane bending coupled with out-of-plane ring deformation 
C-H out-of-plane bending coupled with out-of-plane ring deformation 






